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The one electron spectral functions for the Luttinger model are discussed for 
large but finite systems. The methods presented allow a simple interpretation of the 
results. For finite range interactions interesting nonunivesal spectral features emerge 
for momenta which differ from the Fermi points by the order of the inverse interaction 
range or more. For a simplified model with interactions only within the branches of 
right and left moving electrons analytical expressions for the spectral function are 
presented which allows to perform the thermodynamic limit. As in the general spin- 
less model and the model including spin for which we present mainly numerical results 
the spectral functions do not approach the noninteracting limit for large momenta. 
The implication of our results for recent high resolution photoemission measurements 
on quasi one-dimensional conductors are discussed. 
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I. INTRODUCTION 



The experimental study of quasi one- dimensional conductors can provide a test of the 
pecularities of correlated electrons in one dimension In particular high resolution 

valence photoemission is a very useful tool, as the measured spectra are directly related to 
the one-particle Green's function of the system. Recent experiments of this type give 
a strong indication that the interpretation of spectra requires the inclusion of many body 
effects. As theoretical approaches usually concentrate on the universal behaviour of the 
spectral functions in the extreme low energy regime |3| these experiments provide a stimulus 
to examine also the nonuniversal behaviour of the spectral functions for one- dimensional (ID) 
correlated electrons. 

As first discussed by Tomonaga |S| the problem of ID electrons with a long range in- 
teraction simplifies considerably because it is a good approximation to linearize the energy 
dispersion around the two Fermi points ±kp. In the Luttinger model || an exactly linear 
dispersion is assumed. An exact solution for the Luttinger model was presented by Mattis 
and Lieb 0. The original Tomonaga model and the Luttinger model were compared by 
Gutfreund and Schick || who showed, that the low energy physics in both models is the 
same for long range interaction with a rather weak restriction on the interaction strength. 

The Luttinger model is often studied with the simplification of a zero range interaction 
||. This is sufficient for the discussion of the low energy singularities of the spectra. As the 
interacting ground state of the corresponding model contains holes deep below the Fermi 
level, a direct comparison with a system of nonrelativistic electrons is doubtful, as the 
linearization of the energy dispersion is no longer justified for all relevant energies. We 
therefore study in this paper the spectral properties of the Luttinger model for a finite range 
interaction. 

Even for the spinless case the simplified model with interaction terms only within a 
branch (g^- interaction in the "g-ology" classification ||) the nonuniversal behaviour of the 
spectral function is very nontrivial for /c-values which differ from ±kp by the order of the 
inverse interaction range or more. As the spectral functions of this simplified model can 
be calculated analytically this is probably the most simple nontrivial model of interacting 
electrons for which a complete explicit calculation of spectra can be performed. This solution 
is presented in Sec. [TIT] in two different ways: In a direct approach to calculate the many 
electron eigenstates which enter the Lehmann representation for the spectral functions and 
using the bosonization of the field operators P JT0|JTT|| . In both approaches we first calculate 
the spectra for systems of finite length L and then perform the limit L — > oo. 

For the complete spinless Luttinger model including the ^-interaction terms between 
the branches a recursive numerical method is presented in Sec. |TV] to calculate exact spectra 
for arbitrarily large systems. This procedure gives more insight into the interesting nonuni- 
versal features of the k and u dependence of the spectra than the brute force attempt to 

> 

perform numerically the double Fourier transform of the Green's functions Ga (x, t) which 
themselves have to be calculated involving a numerical integration. For the fc-integrated 
spectral function the direct integration procedure is compared with the asymptotic results 
for very large but finite systems. The approach used for the spinless model is generalized 
to the model including spin in Sec. [V[ Our results and the relevance for photoemission 
spectra are summarized in Sec. |VT[ Spectral moments and the momentum distribution are 
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discussed in Appendices. 

The present paper with its discussion of spectral properties at arbitrary excitation ener- 
gies is complementary to our recent work [TJJ in which results for the universal low energy 
part of the spectra are presented. 



II. SPINLESS LUTTINGER MODEL 

As discussed in the introduction it is necessary to consider the Luttinger model with 
a finite range interaction if one wants the model to work as an approximation to describe 
nonrelativistic electrons not only in the asymptotic low energy region. The original nonrel- 
ativistic electrons are assumed to interact by a two-body potential v (x), i.e. the interaction 
part of the Hamiltonian for a system of finite length L with periodic boundary conditions 
reads 

1 r L r L ,1 

V — - / dx dx'v(x — x')p(x)p(x ) — -v(x = 0)N 

2 Jo Jo 2 

= ^lY.hP-kv{k) + ^N 2 v(k = 0) - ±Nv(x = 0) (1) 

where p(x) is the operator of the electron density pk with k = 2nn/L, n e Z its Fourier 
components and N = f p(x)dx the particle number operator. Here we assume that the two 
body potential v (x) and its Fourier components v(k) are finite at x = respectively k = 0. 
The transition to the interaction term for the Luttinger model occurs by writing pk as 

Pk = Pk,+ + Pk- (2) 

where the pk, a with a = +(— ) are the Fourier components of the operators for the densities 
of right (left) moving particles. The first term on the rhs. of Eq. ([!]) is generalized to 

^ = ^y Yj {94{k)Pk,aP-k,a + g 2 {k)pk,aP-k,-a) (3) 
ZL k^0,a=± 

with the original model corresponding to g2{k) = g&{k) = v(k). The other two terms in 
Eq. (D) are usually neglected. This is justified for the calculation of the spectral functions 
if all frequencies are measured with respect to the chemical potential p. For the explicit 
calculation of p these terms are of importance. One comment on the term involving the 
interaction g±(k) should be made, as it is often dropped. This is only justified for a strictly 
zero range interaction, for which the decomposition used in Eq. ([!]) is not allowed. 
In the spinless model the p q ^ a are given by 

= ^k,oflk+q,a (4) 
k 

where d\ a is the creation operator for a particle of type a and momentum k. With a proper 
normalization 

1 /2 

S = ( x / h,+ for g > (5) 

q -\\q\L) X U,-forg<0 W 
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the density operators obey Bose commutation relations |||6|J7|] 



S, 



1,1 ' 



. 



(6) 



The key to the exact solution of the model lies in the fact, that the kinetic energy can also 
be expressed in terms of the Bose operators M 



T = vp aka\ 



k,a a k,a 



k,a 



v F J2\q\blb q + c(N a ). 

q^O 



(7) 



The additional term c(N a ) involving the particle number operators is irrelevant in the fol- 
lowing and will be dropped. A simple unitary transformation brings the total Hamiltonian 
H = f + V into the form M 



ala 



1 1 



q^O 



where the a q are new boson operators and 



v F (q) = v F ^[l + g,(q)/(27rv F )} 2 - [g 2 (q) / \2ttv f )f . 



(8) 



(9) 



In Eq. (|§D a constant and terms involving particle number operators have been dropped in 
accordance with the discussion of the particle number terms in Eq. ([!]). 

The parts of the one particle Green's function which lead to the photoemission and 
inverse photoemission spectra are 



(10) 
(11) 



These functions can be calculated exactly e.g. by bosonizing the field operators ip a (x,t) 
PJlO|. For finite systems one obtains [|TTJ 



iG l a <] {x,t)e itA = - 



2tt 



ex P tE: 



q>0 



e(-) * eC+) 9 



+2s 2 (g) (cos (gx)e(+) 



ILUnt 



(12) 



where to q = \q\v F (q) and s 2 (q) = sinh 2 (Q q ) with Q q the phase in the unitary transformation 
a q = cosh (Q g )b q — sinh (Q q )P_ q . From (|T2| ) one calculates the relevant spectral function as 



l k,a\ 00 



) 



2tt 



oo poo 



dte 



dxe- ikx e^HG<{x,t) 



(13) 



1 roo roo 

/ dte™ / dxe~ ikx e iflt iG>(x,t) . 

Z7T J —oo J co 



(14) 
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The low energy singularities of the spectral functions have been known for a long time P,|12 

> 

They are obtained by taking the limit L — > oo, in which the Ga {x,t) can be calculated 
analytically in the large x and t limit. The double Fourier transform yields expressions for 
the critical exponents of the threshold singularities of the spectra. 

The aim of our paper is to present results also for the nonuniversal frequency range in 

> 

pcP (k,to). In principle this can be done by direct numerical integrations after performing 

the limit L —>■ oo. As for a finite range potential the frequencies u q have a nontrivial q- 

> 

dependence, the G a (x, t) have to be calculated involving a numerical integration which 
has to be followed by a double numerical Fourier integration. It is numerically difficult 
to obtain the sharp spectral features by this procedure. It also gives little insight into 
the interpretation of the calculated spectra. We therefore take a different approach and 
calculate the spectra for finite systems in such a way that the double Fourier integral can 
be performed analytically. For the simplified model with the ^-interaction only the limit 
L — > oo can be directly read off the results for finite L. As this model already shows very 
interesting nonuniversal behaviour of the spectra we start with a detailed discussion of this 
g 4 -model. 



III. SPECTRAL FUNCTIONS FOR THE SPINLESS G4-MODEL 



As the Hamiltonian for this special model is a sum of Hamiltonians for right and left 
moving electrons it is sufficient to consider e.g. the right moving ones only. In the fermion 
representation the Hamiltonian reads (a n = ak ni + ; p n = Pq n ,+) 



H 4 



2tt 



1 



v F ^2 na\a n + — ^ gi(q n )p nP—n 



(15) 



while the boson representation, apart from an additional term c + (N + ), is given by (b n = b qn ) 



2tt 

H + = T 



v F nblb n + j~Y. n 9i{Qn) (b ] n b n + bJi) 

h>0 47F n>0 



(16) 



In this representation it is obvious that one can read off the exact energy eigenvalues without 
a canonical transformation. The eigenstates are identical to the eigenstates of the noninter- 
acting system. Especially the interacting ground state is given by the Fermi sea | F + (N) ), 
which has the form of a Slater determinant. The general eigenstates have the form 



{m j },N) = H 



1 

m 



1/2 



(if)"' \F + (N)) . 



(17) 



This are linear combinations of electron-hole pair excited states which yield the same value 
of the kinetic energy. It is instructive to write out the states with low excitation energy to 
see the high degeneracy of states with the same kinetic energy. 

The spectral functions for the simplified model can be calculated directly without us- 
ing the bosonization of the field operators. The spectral weights of the delta peaks 
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e.g. in p+(k n , uj) are given by ( {rrij} , N + 1 1 at \ F + (N) ) and can be calculated using 
a+ | F+(N) ) = ala np+l \ F + (N + 1) ), where n F = Lk F (N) / (2rr), and 

/1\ 1/2 

b n= ( - ) E a U«m ( 18 ) 



for n > 1. For n = rip + 1 and = np + 2 the states ajj F+(N) ) are eigenstates: 
at F+1 | F+(JV) ) = | F + (N + 1) ), al F+2 \ F + (N) ) = b\\ F + (N + 1) ). Forn = n F + 1 + n 
with n > 1 the state at | F + (jV) ) — a n F +i+n^n F +i I F-\-{N + 1) ) has overlap to the states 
b[\F + (N + l)), &t&t_ 1 | j p + (iV+ 1)), ^_ 2 |F + (iV + l)), (1/2!) 1 / 2 (6+) 2 S+_ 2 | ^V(AT + 1) > 
etc., i.e. to states | {rrij} ,N + 1) with J2j>iJ m j — n - F° r large n these eigenstates are 
rather complicated linear combinations of electron-hole pair excited states, but the expan- 
sion coefficient of the component a/ nF+1+ii a nF+1 \ F + (N + 1) ) is simple. The square of the 
overlap is given by 

( { mj } ,N + l | ai F+l+h | F + (iV) > I 2 = J} — - = A ^ m ^ ' ^ 

j>i m j- \J ) 

The factor l/rrij\ is due to the corresponding factor in Eq. flT?|) while the factor (l/j) mj 



comes from the prefactor on the rhs. of Eq. (18). This yields for the spectral function 

p > (k nF+l+h} u) = S ^2 j jm j ,n A ({ m j}) S u - H m i w i ( 20 ) 

i m j} 3 V 3 J 

where uij = (2ii/L)(v F + g^qj) / (2n)) . Alternatively this result can be obtained using Eq. 
(|12|) for s 2 (g) = by formally expanding the exponential function and performing the double 



Fourier integral in Eq. ([14]) analytically. We will discuss this procedure in more detail later. 

The calculation of p > (k n ,uj) is therefore reduced to the combinational problem to find 
all decompositions 

mi + 2rri2 + 3rri3 + . . . + nrrin = n (21) 

with rrij e 1N . The solution can be easily produced on a computer, but for large n the 
number of decompositions increases exponentially. For the special g-dependence of g^ used 
in the following 

gA{k)=gS{k 2 c -k 2 ) (22) 

where r c = l/k c is the effective range of the interaction, we will therefore also use a different 
technique. The fc-dependent Fermi velocity v F (k) takes only two different values with this 
assumption 

= (v F = v F (l + gj (2irv F )) for < k < k c 

\v F for k > k c . v ' 

For a finite system this means that g^ik,n) — #4 for 1 < n < n c where n c = Lk c /(2ir) and 
9i{K) = for n > n c . 
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From Eq. (p0|) it follows that p > (k n , uj) is trivial for < k n — k F < k c , as rrij can only be 
different from zero for j < n and all corresponding ujj are given by (2it/Li)vp. This yields 
(k = k — kp) 

p > {k F + k, uj) = 5{uj - vpk) for < k < k c (24) 

i.e. Fermi liquid like behaviour. This is due to the special choice (p^) for g^(k). For k > k c the 
spectral functions are nontrivial. One has to distinguish the intervals mk c < k < (m + l)k c . 
We discuss in the following small values of m and the limit m — > oo. For m = 1 it is still very 
simple to argue in terms of the decompositions in Eq. (|2l"l) . For n c < h < 2n c one can have 
at most one nonzero m/ {mi = 1) for n c < I < h. The remaining "momentum" n — I can be 
decomposed into momenta j which are smaller than n c . Therefore the energy for all these 
decompositions is given by (2tt/L) [Ivp + (n — l)vp\ = (2tc/L) [nv F + (n — l)(v F — v F )]. The 
correponding weight is 1/7. The remaining weight lies in a delta peak at (27r/L)nv F , which 
corresponds to the decompositions of n with nonzero rrij only for j < n c . The limit L — > oo 
can easily be read off and one obtains 

>/, , j- N ©(^ — f F k)Q(vpk + (vp — Vp)(k — k c ) — uS) 



p (k F + k } uj) 



Vpk — uj 



+ 



1 - In (jfe/fcc)] - WFfc) • (25) 



The weight z(k) of the delta peak decreases continously from 1 to 1 — In 2 w 0.307 when 
A; increases from fc c to 2k c . The shape of the spectral weight is shown in Fig. 1(b). If one 
further increases k the calculation of the spectrum using Eq. ( p0|) becomes more and more 
tedious. We will therefore analyse the spectra by another method. 

In the limit k — > oo one might expect that the spectral function p > (kp + k, u) reduces 
to a delta function as for noninteracting electrons. That this is not the case can be seen 
quite generally by calculating the first and second moment of p > (k F + k,u). As shown in 
Appendix A the result for A| = p^i^F + k) — (p^(k F + k)) 2 for the special model of Eq. 
(f22|) for k > 2k c is given in the limit L — > oo by 

A I = l*Z(VF - vp) 2 = \ [p\{k F + k)- vpkf (26) 

i.e. the effective width of the spectrum is independent of k as soon as k is larger than 2k c . 

We now present a method to calculate the spectra using Eq. (|12|). As in the following 
we measure energies with respect to p and momenta with respect to k F we have iG > (x, t) = 
exp [F(x,t)]/L, where for g^ik) given by Eq. 



, N ^ 1 / 2tt . A ^ 1 ( 2tt . A 

F(x,t) = ^ — exp \i—n[x — v F t\ J + ^ — exp I z— n [x — v F t\ J , (27) 

as s 2 (q) = for «2(^) = 0. The second sum on the rhs. of Eq. ( p7| ) is not convergent as it 
stands. There are two obvious methods to overcome this difficulty One can add a factor 
exp (—On) or restrict the sum to n < M, where M(2tt/L) is much larger than the momenta 
one is interested in. Both procedures give the same results for the spectra. Using the first 
method we can write Eq. (|27|) as 
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27T 

F(x, t) = — In ( 1 — exp ( i— [x — v F t + iO] 



+ E — ^exp ^— n [ x ~~ v Ft]J — exp {i~j^ n [ x ~ VFtf^j ■ (28) 

Unfortunatly the finite sums cannot be summed in closed form. We therefore write 

1/L 



iG > (x,t) 



1 — exp {i^Y [x — Vpt + iO 
x 1 1 LL n ex P l ~r m J F - v F t\ 



m=l \j=l 



j\ \ L 



x 



~(l/m)' /.2?r ir 
2^ — ^ — ex P i-j-rnl [x - v F t\ 



, l\ \ L 



= - 2 «^ c) ^p (z^m [x - M] ) E & / ( " c) ex P 'f 1 1* - ) ( 29 ) 

m=0 /=0 

where the affl and &[ nc ^ can be determined iteratively for m > 1, 1 G INq and z = 0, . . . , m — 1 



(m+l) = \ - (m) 

2m+i -| a m(l-j)+i 

j=o J- 



. (m) 



•| m(l-j)+i 
3=0 J' 



The starting values are = \/m\ and 



Using Eq. fl23|) the double Fourier transform can be trivially performed (k F = 0) 

(n c Mn c ) r ( 2tt 



=o 



p>(k n , w ) = £ a^jb^S [u — — [(n — I) v F + lv F ] ) . (32) 



(30) 



# = E("1)Vj! • (31) 

3=0 



Compared to Eq. (|2H| ) this representation avoids the combinatorial problem of Eq. ( |2"T| ) 
and can be used numerically for much larger values of n c than Eq. (p0|). Some important 
information about the a^ c ) and can be obtained analytically. For the case of the 
it is useful to go back to Eq. fl27D , where one can directly read off clq"^ = 1, af a ' = for 



-)•••) Hci Q'i 



(n c ) 



1/i for i — n c + 1, . . . , 2n c + 1 etc.. The limit m — ► oo of a^ c ^ for fixed 



n c follows from Eqs. ([30]) and (|31| ) 

n c / i \ 1/n 

' ( "• , - n ( 

71=1 VC 



^ - n - • ( 33 ) 
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In the large n c limit one therefore obtains — > e~ c /n c « 0.56/n c where C is Euler's 
constant. The behaviour of the for m values of the order of n c follows from Eq. Q27| ) if 
one writes the first term on the rhs. as a difference of a logarithm similar to the first term 
on the rhs. of Eq. (123) and the sum running from n c + 1 to infinity. This yields b\ nc ^ = 1, 



for i = 1, . . . , n c 

t c) = 1 - ( —^—r + + ■■■ + -) (34) 

Vn c + 1 n c + 2 i J v 7 

for n c + 1 < i < 2n c , etc.. 

In this approach it is very simple to obtain the exact analytical result for the spectrum 
in the large momentum limit k^> k c 

n c /-i s 1/m oo y 9 s 

p>(fc n ,cu)^ n - J2t c) t(u--rinv F + l(v F -v F )]) . (35) 
m=l Ve/ *=0 V ^ 7 

We discuss this for the special case n c = 1 and n c — > oo only. For n c = 1 we have 
6| = 1//! and the spectrum is a Poisson distribution with a unit strength parameter. In 
the thermodynamic limit L — > oo, n c — > oo with A; c = n c (2n/L) = const the spectrum is 
given by (5v F = v F — v F ) 

> e ~ c J 1 for < u - WfA; < 5^A; C 

P ( ' W) ~ 5^ X \ 1 - In ) for (J Vj ,fc c < u; - ^A; < 25i; F fc c , ( ° ' ' 

etc.. This behaviour is shown in Fig. 1(d). Fig. 1 summarizes our results for the (74-model 
for v F = 1.2v F , i.e. a repulsive interaction. In Fig. 1(a) all spectral weight lies in a delta 
peak at u = v F k [z = 1). For k > k c there continues to be a delta peak at u = v F k, but its 
weight decreases rapidly with increasing momentum as shown in Figs. 1(b) and 1(c). The 
additional weight lies in the continous part of the spectrum as discussed analytically in Eqs. 
( p5|) and fl36|). In the limit k~^> k c one does not recover the limit of noninteracting electrons. 
The shape of the spectrum becomes independent of the interaction strength as a function of 
the scaled variable (u — v F kp)/(k c 5v F ). A calculation of the spectrum using perturbation 
theory for the selfenergy to low order (e.g. second order) completely fails to give the correct 
shape of the spectrum. 

Another quantity of interest is the total spectral density per unit length 

/•OO rjh 

p>{u)= / - P >(k,u) (37) 

1 r°° 

/ d tG>(0,t). 

Z7T J—oo 

For the special ^-interaction ( p2[ ) the function p > (k,u) has been given analytically in Eqs. 
(^4j) and (|25|). The /c-integration can be simply performed and yields 

(1 for < uj < v F k c 

1 + In (-^) for v F k c < 00 < v F k c (38) 
1 1 ford F k c <u< 2v F k c , 
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etc.. Alternatively the function F(0, t) in Eq. ([27]) can be calculated numerically in the limit 
L — > oo and G > (0, t) is Fourier transformed numerically. The results in both approaches 
agree and are shown in Fig. 2. The latter method can be performed for arbitrary interactions 
g±{k). Results for g±{k) = #4 exp (— \k\/k c ) are also shown in Fig. 2 for a positive and a 
negative value of g±. For the case of a repulsive interaction there is a depletion of the total 
spectral weight at low frequencies. Due to a sum rule discussed in Jl3| the missing weight has 
to show up in another frequency range. For a model without an upper cut-off in momentum 
space the missing weight is pushed to infinity. If a cut-off is included the missing weight 
appears at the upper end of the spectrum, as the spectral weight in Fig. 1(d) has to be 
compared to a delta peak at vpk max for noninteracting electrons. 

In this section the simple (74-model has been discussed quite at length, as the generaliza- 
tion to the complete model including #2 terms has to be performed mainly numerically. 



IV. SPECTRAL FUNCTIONS FOR THE GENERAL SPINLESS MODEL 



The method to calculate G > (x, t) presented in Eq. ( p9|) is generalized to the full spinless 
model introduced in Sec. [TI[ We specialize to interactions of the type described in (p2|), i.e. 
we also assume g2{k) = g2&{k 2 — k 2 ). Then using Eq. fll2"D the expression for F(x, t) in 
is generalized to (s 2 = s 2 (q n ) for 1 < n < n c ) 



F(s,f) = E 



n=l 



n 



. 1 2tt A 2 / 27T 

1 + s ) exp ( 2— ^ [x — Vpt\ J + s exp I —i—n [x + vpt\ 



2tt 



2s — exp ( i—n [x — vpt] 



In 



( 2tx . A 
1 — exp y~^ n F — v Ft + *0J I 



(39) 



Again iG > (x,t) = exp [F(x,t)]/L is written as a product of power series 

1 / 00 / 2n \ \ 

iG > (x, t) = — I a^ 1 exp ( i—m [x — vpt] ) ] 

L \m=0 ^ L 'J 

1 



x i^2b\ n ^ exp (i-yl [x - v F t] 
\i=o 

( °° / 27T \ \ / ™ c 1 ^ 

x ( 4" c) exp ( -i—r [x + ) ) exp -2s 2 ^ - 



\r=0 



n=l 



n 



(40) 



where the expansion coefficients are determined iteratively as in Sec. III]. For m > 1, I G N 
and z = 0, . . . , m — 1 one obtains 



(m+l) 



' (-l/ m )i (m) 

im+i -| a m(l-j)+i 

3=0 J' 

1 ([l + s 2 }/my u(m) 

m(l-j)+i 

3=0 



a 



, (m+l) _ \ - 



(41) 



' 's 2 jmf 



-lm+i ■ f u m(/-j')+i • 

j=0 J- 
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The starting values are b$ = (1 + s 2 ) m /m\, c$ = s 2m /m\ and is given by Eq. (^TJ). 
Using Eq. ( f4"0| ) the double Fourier transform can be simply performed (A = exp (X)«Li 



oo n+r 



27T 

x () ( jj — — [(n + r — j) + (r + j) f f] 



(42) 



If we write (n + r — j)f_F + (r + j)?Jf = nt>F + j(vp — vp) + r(vp — Vp) it is obvious that 
p > (k n , uj) = for uj < k n vp is guaranteed only for vp — vp > 0, i.e. for repulsive interactions 
if 9a = 92- 

The coefficients affi^ are the same as in Sec. [TTI[ The behaviour of the b]^ c > and for 
m < n c follow from the identity 



ex P ( 7 £ 



71=1 



n 



[1 - z) 7 exp -7 E 



n=n c +l 



77 



oo / m 



i+E n 

m=l W=l 



1 + 



7 



OO 

x exp I —7 E — ^ 



n=n c +l 



77 



(43) 



For 1 < m < n c the expansion coefficients are therefore given by the first factor on the rhs. 
of Eq. PD, i.e. 



e» = ni+- , 

1=1 v 3 j 

s 2 -l 



— > const X 777 



(44) 



,(ra c ) 



n 1+ 

3=1 \ 



n c >m>l 



J 



const x 777 



s a -l 



This power law behaviour of the coefficients and is responsible for the low energy 
power law singularities of the spectral functions and the power law behaviour of the momen- 
tum distribution n(k) discussed in Appendix B. To demonstrate this we consider momenta 
1 <^ n <^ n c and frequencies of the order vp(2iin/L) in Eq. ([Hp. As discussed before Eq. 
(|33|), = 5 m fi for m < n c . Therefore the spectral function simplifies in this regime to 



p>(k n , u ) = E c^b^lS L-^[n + 2r] v F ) 

r>0 \ L / 



(45) 



i.e. the spectrum consists of delta peaks at u — vpk n + Vp(Aitr/L). For 1 < r + n C n c we 
can use the asymptotic form of the coefficients in Eq. ( f44"l) at r = (uj — vpk n ) / (Anvp/ L) and 
77 + r = {uj + v pk n ) / (Anil p / L) to obtain the weights of the peaks. In the limit L — > 00 this 
yields 



p > (k n , uj) ~ 9(c<j — vpk){uj — vpk) s 1 (u> + t^fe)* 



(46) 
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i.e. the well known asymptotic behaviour for k <^ik c and ui — vpk <C vpk c [j4|Jll~||. 

In the opposite limit k 3> k c the coefficients a^l-j in Eq. fl42|) can be replaced by the 
constant introduced in Eq. fl33f ). Near the threshold at u — vpk n Eq. (fHf) can be used 
again for the and b™ . Performing the integrations in the limit L — > oo the finite step 
at threshold in Eq. fl36| ) is replaced by a power law behaviour (u — vpk) 2s2 . For arbitrary 
values of fc n and u the spectral function p > (k n ,uj) has to be calculated numerically. Results 
for the same k values as used for the ^4-model in Sec. | IT] and gi (k) = g^{k) are shown in 
Figs. 3(a)-(d). The figures shows that the delta peak at to = vpk is changed into a power law 
behaviour and one could analytically show that the critical exponent for u \ vpk is given 
by s 2 — 1 as in Eq. fl46|). The shape of the continous part of the spectra of the (74-model is 
modified quite considerably for the parameters used (s 2 = 1/8 corresponding to vf = 2vf)- 
In contrast to the (74-model p > (k,u) 7^ for k < and small \k\ [ llj, but the spectral weight 
for negative k decreases rapidly for increasing \k\. 

In Fig. 4 we show the integrated spectral weight for g±{k) = g2(k) for the step model 
(HI) as well as the exponential model g&(k) = g2{k) = gexp(—\k\/k c ) used also in Fig. 2. 
At low frequencies the power law behaviour proportional to u 2s leads to a suppression of 



spectral weight [14]. For the case of an attractive interaction this leads to a peak in p > {uj). 



V. GENERALIZATIONS FOR THE MODEL INCLUDING SPIN 

For the Luttinger model including spin the discussion in Sec. |TI| applies up to Eq. ([|). If 
the density operators p q ^ a are decomposed into a sum of particle-hole operators an additional 
spin summation occurs. It is then useful to define charge and spin operators 



Pq,a — Pg,o,T + Pq,a,l {^7) 
®q,a = Pq,a,1 Pq,a,i • 

With a normalization which differs by a factor of \[2 the analogous definition to Eq. ([5]) 

reads 

q ' C \\q\L) \p 9 ,-forg<0 



q\L I 1 (t q - for q < 



These operators describe independent boson degrees of freedom. Again the kinetic energy 
can be expressed in terms of the boson operators and particle number operators 

f = v F £ |g| (b\Xc + KAs) + C (N) . (49) 

q^O 

For the spin independent interaction ([[]) the spin degrees are not renormalized by including 
the interaction, i.e. VF, s {q) = vf- For the charge degrees of freedom the problem to find the 
exact eigenstates is equivalent to the spin-independent problem. The additional factor \/2 
in Eq. © modifies Eq. (|) to 
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VfM = v F \j[l + g±{q)/{nv F )f - [^ 2 (g)/(vrt;/)] 2 . (50) 



For the propagators the changes are more dramatic. If one denotes the propagator for the 

> 

spinless model by Ga {x,t; g 2 , g±) one obtains using the bosonization of the fermion field 
operators 

iO$ (x, t; g 2 , g 4 ) = [iG { I <] (x, t; 0, 0>'G# (x, t; 2g 2 , 2gSj ■ (51) 

For the ^4-model the more intuitive approach to calculate the spectral functions via the 
Lehmann representation described for the spinless model in Sec. |T| can be generalized and 
leads to a combinatorical problem slightly more complicated than in Eq. (|21|) . 

In the following we restrict the discussion to the step model (|22|) . For the (74-model the 
spectral function has the same form as (B2j) 



p>(k n , w ) = £ at\A nc) S (« - -f [nt* + / (^ - 1*)] . (52) 

The coefficients /^""^ are obtained from the power series with the coefficients b^(2g±) by 
taking the square root, i.e. 

A K) = (V c) - £ /5ft } /5] nc) j / c) ) • (53) 

The coefficients a^j are obtained from = Yl™=o a n"^ (2^4) by the same procedure. For 
the (?4-model including spin the spectral density differs from a simple delta peak already for 
k < k c . For I < n c the coefficients &i and can be given explicitely as = 1 = 6^ c ^ 
for m < n c . For 1 < I <C n c one has 5j = /^""^ ~ Z -1 / 2 . In the limit L — > 00 this yields 
for < < fc c and -up = vp tC > vp 

p > (k, to) = const x ([a; — Vpk] [vpk — c^]) -1 / 2 0(co> — vpk)Q(vpk — (54) 

as discussed previously []TT|, IB , TT ] . Note that this non- Fermi liquid like behaviour occurs, 
but one still has n(k) = for k > kp = 0. The behaviour of the spectral function is shown in 
Fig. 5 (dotted curves). Even for k ^> k c there remains a square root singularity at threshold 
which is "intermediate" between the delta peak expected from the first factor in Eq. (|5T|) 
and the plateau from the second factor. 

For the general model including spin the spectral function follows from Eq. fl5"T|) using 
Eq. © 

00 n+r 

p>(k n , u) = A~ s2 £ £ a^l^fH na) 

r=0 j=0 

x5(uj-^-[(n + r-j)vF + (r + j) v F ]j (55) 

where the coefficients /5>j and 7^ nc - ) follow from the coefficients b^ c \2g2,2gi) and 
(2# 2 , 2g 4 ) defined in Eq. ( fiTf ) by the procedure ( |53D to take the square root of a power 
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series, while the cq follow from the = YhLi a i (2<?2j 2^4) in the analogous way. In 
the low energy regime Eq. ( pH ) leads to the power law behaviour discussed in detail in 
references |1]J and JET]. This and the spectral behaviour for larger values of k is shown in 
Fig. 5 (solid curves), where the results of the full model are compared to the ^4-model. For 
the value of s 2 used the nouniversal features of the spectra are rather similiar, i.e. they can 
be largely understood examining the much simpler g 4 -model. 

Results for the integrated spectral density p > {yj) are qualitatively the same as the cor- 
responding curves for the spinless case shown in Fig. 4. 



VI. SUMMARY 

In the preceeding sections we have presented results for the spectral functions p^(k,uj) 
and p%{uj) relevant for inverse photoemission. The corresponding functions to describe 
photoemission are given by the mirror images of the curves presented as p%(kp + k,u) = 
P+(k F — k, —uj). In the Luttinger model with its linear energy dispersion the value of the 
Fermi momentum is irrelevant and has been put to zero in Sec. ITT|-|V|. If we want to 
describe nonrelativistic electrons this is no longer the case as kp is proportional to the 
electron density. In order to apply results using the Luttinger model one has to keep in 
mind that the linearization procedure is only allowed for sufficiently long range interactions, 
i.e. k c <C kp. 

In Ref. |J angular integrated high resolution photoemission data of quasi one-dimensional 
conductors are presented which show a depletion near the Fermi energy and a rather broad 
peak about one eV below the Fermi level. As these spectra are related to p < {oj) the peak 
below the Fermi level reminds one of the peak in one of the integrated spectral functions in 
Fig. 4. But this spectrum corresponds to an attractive interaction. Such a peak was first 
discussed by Suzumura who obtained it for the case g 2 > 0, g 4 = 0, i.e. a repulsive g 2 - 



interaction but neglecting the ^-interaction. As one can see from Eq. (^) or ( |5DD this leads 
to Vp(k — ► 0) < vf, i.e. an increase of spectral weight with decreasing |a>| until the power 
law factor u a (a = 2s 2 for the spinless model, a = s 2 for the model including spin with spin 
independent interaction) leads to suppression and the peak emerges. As discussed in Sec. |I| 
the physical model corresponds to g 4 = g 2 , i.e. it is unphysical to neglect the (^-interaction 
for finite range interactions. Therefore the experimental peak below the Fermi level cannot 
be explained as a Luttinger liquid feature for repulsive interactions. In order to explain the 
experimental depletion near the Fermi level a suprisingly large value of the exponent has to 
be assumed |L4| . 



In this paper we have presented the first detailed study of the nonuniversal spectral 
properties of the Luttinger model and have shown both analytically and numerically that a 
suprisingly rich variety of spectral features can emerge. 
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APPENDIX A 



In this appendix we calculate the moments 



fin(k) = / uj n p (k, u)duj 



'd n G>(x,ty 
dT- 



~ tkx dx 



(Al) 



t=o 



of the spectral function p > (k,u) for n — 0,1 and 2. As iG > (x,t) = exp [F(x,t)]/L we have 
to calculate the derivatives of F(x, t) 



G > (x,0) = F(x,0)G > (x,0) 

G>(x, 0) = (P(x, 0) + F 2 {x, 0)) G>(x, 0) . 



(A2) 



The calculations are straightforward only for the ^4-model. In the spinless case the function 
F(x, t) takes the form 



F(x, t) = — ex P (i~r n l x ~ v F( n )t] 

n>l 71 ^ L 



(A3) 



Therefore F(x, 0) = In (1 — exp 1^ (x + iO) ) , i.e. G > (x, 0) is identical to the noninteract- 
ing Green's function 



1 °° / 27T \ 

iG > (x, 0) = — ^2 ex P ( i-j-n [x + iO] I 
L n =o ' ^ L ' 



(A4) 



The spatial Fourier transform in (|A1|) can easily be performed using (|X2-X4) and yields 

(v F (i) = v F + g^ih)/ (27r)) 



Pi (k n ) 



L m=l 

vUkn) ~ (Pi(k n )) 2 

n s 2 1 n n 

7 EE (v F (i)-v F (j)f ■ 

L> ' Z j= l l= n-j 



(A5) 



(A6) 



For the step model 



this leads for n > 2n r to 



A>= ^ 



2tt\ n c (n c + 1) 



(v F - v F f . 



(A7) 



In the limit L — > 00 this reduces to Eq. (|26|) . 

For the general spinless model already the zeroth moment p^ (k n ) = 1 — n(k n ) is non- 
trivial. It is discussed in Appendix B. As the expressions for the first and second moment 
are rather lengthy we do not present them here. 

For the model including spin we restrict ourselves to the ^4-model. If we again write 
iG>(x,t) = exp [F(x,t))/L the function F(x,t) is given by 
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F(x,t) = - *y — ^exp (^~j~ n \ x ~ VF(n)t]j + exp {^—j~ n [ x ~ VFtlj^J ■ (A8) 



Therefore G^(x, 0) equals the expression on the rhs. of Eq. (|A1|) and the spatial Fourier 
transform can be performed as in the spinless case. The first and the second moment are 
given by 



2n n 

Vl.A k n) = jE {VF{™>) + V F ) /2 
^ m=l 



(A9) 



q y 2 i n n 

^ j i j=1 i=n _j 

(n y 2 n 
T ) Ej(Mj)-^) 2 /4. 



(A10) 



As now = vp + giik^/n the first term on the rhs. of ( |A10| ) is identical to the result 

for the spinless model. For the step model the additional term is the only contribution to 
the width of the spectrum in the interval < k < k c . 



APPENDIX B 



In this Appendix we present a short discussion of the momentum distribution in the 
ground state for finite systems and in the limit L — > oo. For the special interaction g2{k) = 
g^{k) = gQ(k 2 — k 2 ) the Green's function G > (x, 0) for the spinless model can be written as 

iG>(x,0) = ± ^(^f[N + 1] x) 



Ll-exp fe[z + i0] 



x exp s' 



y - 



■/ L)nx 



+ e 



-i(2ir I L)nx 



n=l 



A~ 2s2 exp (ikp [N + 1] x) 
L 1 - exp [x + i0] 



y exp ii—mx 

Vm=0 ^ L 



2vr 



x ^2^c ; exp y—i—lx 
with the coefficients given in Eq. (41) and 



A = exp ( y 



n c >l C 



e n r 



Vn=l 



n 



(Bl) 



(B2) 



The momentum distribution n(k) follows from the spatial Fourier transform of iG > (x, 0) 

(B3) 



1 " n f [n F + 1 + n] )= y y dfrM^S, 

V ^ J m=0 1=0 



l+ii,m 
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where = YaLo c i ■ The asymptotic behaviour of d^ can be obtained from a compar- 
ison of the power series in ( |B1[ ) for x = 

oo 

d (nc)---^ c (nc) =A a» (B4) 
1=0 

Using (|B4] ) the large momentum behaviour n ( ^ + 1 + n]J — >■ for n — > oo can be read 
off Eq. (P3|). As 1/2 — n(A;) is symmetric with respect to (2ir / L)(n F + 1/2) it is sufficient 
to consider values h > 0. 

In order to obtain the power law behaviour of the momentum distribution for 1 < n < n c 
it is useful to calculate the finite differences A s = n [n F + hfj —n [n F + 1 + hfj . They 
are given by 

oo 

Aft = A~ 2s2 c\ nc) c\li . (B5) 

1=0 

The singular contribution can already be obtained by restricting the summation in (|B5| ) to 
values n + l < n c . This simplifies the discussion as an analytical expression for the c 2 ( f44[) 
is available. If one uses the asymptotic form ~ m s _1 , Eq. ([B"5|) involves a summation 
over (1/ [/ (/ + n)]) 1_,s which in the thermodynamic limit for 2s 2 < 1 contains a singular 
contribution proportional to (l/h) l ~ 2s . We therefore recover the well known result j?],|8| 



n(k) - 1/2 ~ sign(A; F - k) \k - k F \ 2s2 . (B6) 



If one wants to determine the exponent of dn/dk numerically from Eq. ( |B"5| ) one has to go 
to very large values of n c . 



17 



REFERENCES 



[1] D. Jerome and H. J. Schulz, Adv. Phys. 31, 299, (1982) 
[2] J. Solyom, Adv. Phys. 28, 201, (1979) 

[3] B. Dardel, D. Malterre, M. Grioni, P. Waibel, Y. Baer, C. Schlenker and Y. Petroff, 

Europhys. Lett. 19, 525 (1992) 
[4] A. Luther and I. Peschel, Phys. Rev. B9, 2911, (1974) 
[5] S. Tomonaga, Prog. Theo. Phys. 5, 544, (1950) 
[6] J.M. Luttinger, J. Math. Phys. 4, 1154, (1963) 
[7] D.C. Mattis and E.H. Lieb, J. Math. Phys. 6, 304, (1965) 
[8] H. Gutfreund and M. Schick, Phys. Rev. 168, 418, (1968) 
[9] see Ref. 2 for a review 
[10] F.D. Haldane, J. Phys. C14, 2585, (1981) 

[11] V. Meden and K. Schdnhammer, Phys. Rev. B46, 15753, (1992) 
[12] A. Theumann, J. Math. Phys. 8, 2460, (1967) 
[13] Y. Suzumura, Prog. Theo. Phys. 63, 51, (1980) 

[14] K. Schonhammer and V. Meden, J. Electron Spectroscopy and rel. Phenomena (in press) 
[15] A. Luther and V.J. Emery, Phys. Rev. Lett. 33, 2460, (1967) 
[16] H.C. Fogedby, J. Phys. C9, 3757, (1976) 
[17] J. Voit, Phys. Rev. B, to be published 



18 



FIGURES 



FIG. 1. Spectral function p > (k,u) of the spinless (74-model as a function of uj/{vpk c ) for 
vp = I.Ivf and different momenta as indicated in the figures. The arrows represents delta peaks 
with weight z(k). 

FIG. 2. Total spectral density p > (uj) of the spinless (74-model as a function of uj/(vpk c ). The 
solid curve shows the result for the step model with vp = \.2vf, the dotted curve the result for the 
exponential model with g4/(27rvp) = 0.2 (repulsive interaction) and the dashed curve the result 
for the exponential model with g4/(27rvp) = —0.2 (atractive interaction). 

FIG. 3. The same as in Fig. 1 but for the full #2 = 9 a spinless Luttinger model for vp = 2vp. 
The small oscillations in the full curves are a finite size effect. The dotted curves shows the 
continous part of the related spectral function for the spinless (74-model. 

FIG. 4. The same as in Fig. 2 but for the full §2 = 54 spinless Luttinger model. The solid curve 
shows the result for the step model with vp = 2vp, the dotted curve the result for the exponential 
model with gj (ttvp) = 3 (repulsive interaction) and the dashed curve the result for the exponential 
model with g/{-Kvp) = —0.49 (atractive interaction). 

FIG. 5. The same as in Fig. 1 but for the model including spin. The dotted curves present the 
results for the (74-model including spin and the full curves the related results for the full (72 = 9a 
Luttinger model including spin. The small oscillations are again a finite size effect. 



19 



